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Abstract 

We consider a single spin in a constant magnetic field or an anisotropy field. 
We show that additional external time-periodic fields with zero mean may 
generate nonzero time-averaged spin components which are vanishing for the 
time-averaged Hamiltonian. The reason is a lowering of the dynamical sym- 
metry of the system. A harmonic signal with proper orientation is enough to 
display the effect. We analyze the problem both with and without dissipa- 
tion, both for quantum spins (s = 1/2, 1) and classical spins. The results are 
of importance for controlling the system's state using high or low frequency 
fields and for using new resonance techniques which probe internal system 
parameters, to name a few. 
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Usually nonzero averages of observables, which would be expected to be zero by symmetry 
considerations, are generated either by constant external fields, or by internal interactions 
which may lead to phase transitions. However as we will show below such a situation 
is also possible if we use time-periodic fields with zero mean. The general idea behind 
the following results is purely symmetry related, and thus it seems to be worthwhile to 
understand the mechanisms which may lead to nonzero averages if such fields are applied. 
This work is motivated by a recent paper M where similar ideas have been used to explain 
the phenomenon of directed currents in driven systems. The essence of the present paper 
is that we can lower the symmetry of a given dynamical system by applying time-periodic 
fields with zero mean, i.e. that the time-averaged Hamiltonian displays symmetries which 
would imply zero averages for corresponding observables. It will be the symmetry breaking 
in the temporal evolution which induces nonzero averages. 

Let us start our considerations with a model describing an s = 1/2 spin in a constant 
field h z = 2 directed along the z-direction and a time-periodic field 2h x (t) with period T and 
zero mean directed along the x-direction. The Hamiltonian is given by H = h z S z + 2h x (t)S x 
(here S x , y ,z are the spin component operators related to the corresponding Pauli matrices, 
e.g. [0]). For the moment we assume that <C 2 and the frequency to = 2n/T <C 2. In 

that case we can use the adiabatic approximation and neglect Zener transitions. The two 



eigenvalues of H for a given value of h x are \± = ±y 1 + h x . The expectation value for S x 
in these states is given by 

(s*) = hx ■ (1) 



Now we assume that the spin is in any of the two states. Slow variation of h x in time 
will keep the system in that state. Let us average (S x ) over one period of oscillation. 
Because (S x ) is odd in h x , we will obtain nonzero time averages for the x-component of 
the spin if e.g. Jq h x dt ^ 0. This is possible if h x (t) contains several harmonics (SH), e.g. 
h x (t) = hi cos(ut) + h 2 cos(2ut + £) (see also Q). In that case in lowest order in hi, h 2 we 
obtain (S x ) = — ^h\h 2 cos£. We conclude this example with stating that it is possible to 
generate a nonzero average S x spin component by applying a permanent field in z-direction 
and a time-periodic field with SH and zero average in x-direction. 

Let us relate the results from the example given above to symmetry considerations. The 
Hamiltonian H should be a periodic function of time H(t) = H(t + T). Instead of solving 
the time-dependent Schrodinger equation, which would bring us to the analysis of unitary 
Floquet matrices [[§ , we follow the density matrix approach, which is suitable since we want 
to average over different initial conditions and are thus facing the dynamics of mixed states. 
The density matrix p satisfies the quantum Liouville equation 

^=i[H,p)-u(p- Pp ) (2) 

where [A, B] = AB — BA, pp is some equilibrium density matrix parametrized by the inverse 
temperature (3 and v is a phenomenological parameter measuring the coupling strength of 
the system described by H to some environment. Note that v is the characteristic inverse 
relaxation time of H in the environmental bath. 
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Let us further define H = 1/T J T H(t)dt and H^t) = H(t)-H . Note that H^dt = 
0. Then we may choose pp = ^e~^ H ° with Z = Tr(e _/3H °). We define the value A(t) of an 
observable characterized by the operator A as A(t) = Tr(Ap(t)). The time average of 
A(t) shall be defined as A = lim^^oo h Jq A(t)dt. The averaged attenuation power (the 
rate of energy transfer from the time-periodic field to the heat bath) is given by W = 
u(H -Tr (H oP p)). 

We chose the relaxation term in (Q) in an oversimplified form. There are many theories 
which exploit different concrete relaxation mechanisms (e.g. [|J and references therein). The 
reason for choosing (|2|) instead is that it allows to discuss the following symmetry breaking 
without entering the details of the concrete dissipation mechanism. In other words, we 
deliberately choose the simplest dissipation term which conserves all symmetries of our 
dynamical system except time reversal. 

Equation (0) is a linear equation for the matrix coefficients of p with inhomogeneous 
terms due to pp. The general solution is given by a sum of the general solution of the 
homogeneous equation (put pp = in @) and a particular solution of the full equation. 
Since the homogeneous solution for v = is given by some unitary time evolution, v > 
will cause all solutions of the homogeneous equation to decay to zero for infinite time. For 
t ^> 1/v any particular solution of the inhomogeneous equation trends to a unique time- 
periodic solution - the attractor of (0). This allows us to choose any (reasonable) initial 
condition p(t = 0). If H, p(t = 0) and pp are invariant under certain unitary transformations, 
it immediately follows that p(t) keeps those symmetries, and consequently the attractor will 
have the same symmetries too. For large temperatures pp is approaching the unity matrix 
(up to some factor). Consequently in that limit, whatever the time dependence of H(t), the 
solution of (H) will approach pp. Finally we note that due to Trp^ = 1 any choice of p{t = 0) 
with Trp(t = 0) = 1 implies Trp(i) = 1 for all t. 

Let us consider (0) for 

H = h S z + h(t)(aS x + 1 S z ) (3) 

where a = sin(0) and 7 = cos(0). This model describes a spin in a constant magnetic 
field pointing in the ^-direction, under the influence of an additional time-periodic field 
h(t) = h(t + T). This oscillating field should have zero mean: J h(t)dt = 0. Let us define 
h(t) having T a symmetry if h(t) = —h(—t) = h a (t), T s symmetry if h(t) = h(—t) = h s (t), 
and T s h symmetry if h(t) = —h(t + T/2) = h s h(t) (note that in the two first cases any 
argument shift is allowed, so that e.g. h(t) = cos(t + p) posesses all three symmetries). For 
a monochromatic field (MCF) h(t) and <ft = tt/2 (|3|) is the classical setup for performing 
magnetic resonance (MR) experiments 0,0. 

For the s = | case the spin component operators are given by the Pauli matrices: 
Sx, y ,z = \&x, y ,z- The density matrix p has three independent real variables. Using the 
variables S x ^ y ^ z we find 

S x = (h +>yh(t))S v -vS w (4) 
S y = ah(t)S z - (h + qfh(t))S x - uS y (5) 
S z = -ah(t)S y - u(S z - C) (6) 

where C = l/2tanh(/i /5/2). Note that the obtained set of equations for v = is equivalent 
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to the Heisenberg equations for the operators S XiVjZ and thus also to the equations of motion 
for a classical spin. In fact is a particular case of the Bloch equations 0,0. 

Let us discuss the symmetries of (f|)-(|6|) which conserve H , i.e. S z — > S z . Consider 
the case 7 = 0: if h{t) = h sh {t) then a symmetry operation Qi is S x — > —S x , S y — > 
—S y , S z — > S z , t — > t + T/2 . If Q\ holds we conclude that S x = S y = 0, while S z 
may be nonzero. Consider 7 = and v — 0: if h(t) = h a (t) then a symmetry operation 
Q2 is S x — > —S x , S y — ► S y , — » <S Z , t — > — £ . If Q2 holds it follows = 0, 
while S 1 ^ may be nonzero. Finally for z/ = and h(t) = h s (t) a symmetry operation Q 3 is 
S x ^ S x , S y — > — 5^ , 5* z — > 5 Z , t — >• — t . If Q3 holds it follows = 0, while S 1 ^ may be 
nonzero. 

Let us note some consequences. If we choose h(t) = hicos(out), then the classical MR 
setup with 7 = (Qi) yields nonzero values for S z only p. If the probing field is not 
perpendicular to the z-axis (7 7^ 0), nonzero values for S x and S y appear as well. S y will 
vanish in the limit of zero coupling to the environment v — » (Q3), so that this average 
can be used to measure the coupling strength. Applying e.g. h(t) = hi sm(out) + h 2 sin(2u;t) 
(having h a symmetry but not h s h and h s one) we can suppress the value of S x relatively to 
S y for 7 — > and v — > keeping S y finite (Q2)! 

Analytical solutions to (H)-(H) can be found e.g. for large v ^> 1. Expanding in 1/v and 
averaging over time we find in lowest orders 

S x = C ai {h 2 )^- - Ca(-j(hh) + ^h 2 (h 2 ) + 

(a 2 + 37>o</> 3 > + 7(7 2 + « 2 )(^ 4 )) \ + O(-) (7) 



z/ 4 V 5 ' 



S y = —Col 



2 1 h (h 2 } + ( 1 2 + a 2 )(h 3 )}±- + O(\) (8) 



z/ 3 z/ 5 ' 



where (f(t)} = ^ f$ f(t)dt. It is easy to cross check that all symmetry statements from 
above are correct. Nonzero values for (h 3 ) can be obtained e.g. with h(t) = hi sin(a>£) + 
h 2 sm(2ut + £) for £ 0, tt (see also ||). 

In Fig.l we show the dependence of on for = V2 COSLVt, (f) = 7r/4, /z — 3, 

z/ = 0.1 and (3 = 10. The time-periodic field has a large amplitude compared to typical 
MR setups 0. This causes the S z curve to show a rather broad peak at uo pa h - the 
position of the expected MR resonance. However we also observe sattelite peaks at lower 
frequencies which are clearly related to the variations of nonzero S x , y (for convenience these 
averages are scaled by a factor of 10 in Fig.l). In fact the positions of the sattelite peaks 
are subharmonics of the main resonance. The dependence of S x and S y on u shows rather 
complex structures. We find that typically the dependence of these averages on u becomes 
oscillatory for small u ho, whereas large uj values yield smooth decay curves. Note also 
that these averages stay nonzero down to small frequencies in accord with the adiabatic 
example from above. Also important is to notice that the fluctuations of S x and S y around 
their mean values may happen with amplitudes being one order of magnitude larger than 
the mean values (see inset in Fig.l ). 

The above results hold also for larger spins. To show that they also hold for internal 
anisotropy fields rather than external fields, we consider a spin with s = 1 and the Hamil- 
tonian 
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H = S% + h{t)(aS x + yS g ) 



(9) 



which describes a spin with an anisotropy along the z-axis (S%) under the influence of an 
external magnetic field h(t) parallel to the xz plane. The magnetic field is again time- 
periodic with period T and has zero mean. The 3x3 hermitian density matrix p has 8 
independent real parameters. Since H in (H) is a real symmetric matrix, we can define 
p = R + il where R is a real symmetric matrix and I a real antisymmetric one. Noting that 
also pp is a real diagonal matrix, (||) can be rewritten as 



dR 
~dt 



[H, 1} - v{R - pp) 
dl 



(10) 
(11) 



It follows S x = V2(R(1,2) + R(2,3)), S y = ~V2(I(1, 2) +1(2, 3)) and S z = R(l, 1) -R{3, 3) 
||. Using the abbrevations P x = y/2(R(l,2) - R{2,3)), P y = y/2(I(l,2) - 1(2,3)), P z = 
R(l, 1) + R(3, 3), R 13 = V2R{1, 3), J 13 = \/2/(l, 3), R 22 = V2R(2, 2), D' 1 = 1 + 2e~P and 
= 2 + e' 3 the equations of motion become 



S x = -P y + -yhSy - uS x 
P x = Sy - yhPy + \/2ahIi 3 - vP x 
S y = —P x — yhS x + ahS z — uS y 



P y = S x + yhP x + ah \V2R 22 -P z - V2R 



H3 



Rl3 — 
R22 = 



S z = ahSy — vS z 
P z = ahPy -v[P x - 2F) 
-27M13 + V2ahP y - uR 13 
= 2yhR 13 - V2ahP x - ul 13 
\f2ahPy - v \R 22 - V2D 



(12) 



These equations conserve the trace Trp = P z + R 22 j\[2 = 1. 

Now we can discuss the symmetries of (O) which change the sign of S. Two of them 
hold only for v = 0. First, if h(t) = h a (t), then the equations are invariant under change 
of sign of the variables t, S x , S y , S z (leaving all other variables unchanged). A second case 
takes place if h(t) = h s (t). Then changing the sign of t, S y , P y , I\ 3 (leaving all other variables 
unchanged) is an operation which keeps equations (0) invariant. These two cases imply 
that if h(t) is antisymmetric, then for vanishing dissipation v — > S Xjy>z — > 0, while for 
symmetric h(t) the same limit provides a vanishing of the y-component only S y — > 0. 

For the general case v 7^ two more symmetries may take place. If 7 = (the field 
h(t) acts perpendicularly to the anisotropy axis z), changing the sign of S y , S Z ,P X , I\ 3 (and 
keeping all others) leaves (|12j) invariant. Finally if h(t) = h s h(t), the shift t — > t + T/2 and 
simultaneous change of sign of the variables S x , S z , P y , I 13 do not change the equations. It 
follows that S y = S z = for 7 = and S x = S z = for h(t) having shift symmetry. 

It is interesting to note that for a MCF h(t) = cosut and v 7^ 0, 7 7^ the spin will point 
on average in y direction, i.e. perpendicular to the plane spanned by the driving field and 
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7 = 0: 
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the local anisotropy axis! In Fig. 2 we plot the dependence of S y on u> for this case ((3 = 10, 
v = 0.1, 7 = a = 1), which confirms the symmetry considerations. Note that S x and S z are 
less than 10~ 8 as found in the numerical studies. 

To conclude this case we remark that it is again an easy task to perform expansions in 
\jv for large v values as shown above for the s = 1/2 case. The resulting expressions also 
confirm the symmetry considerations. 

So far we have discussed the results for quantum spin systems. It is also possible to 
analyze corresponding classical systems. E.g. the classical equations for (|j) are given by 

s x = -2s z Sy + jhs y (13) 
s y = 2s z s x + h (as z - 75^) (14) 
s z = —ahs y (15) 

Let us discuss the symmetry properties of (|I~3|)-(|T5|). We denote on the left part the condition 
and on the right part the symmetry operations which leave the equations of motion invariant 
(note that we list only those variables which have to be changed): 

, s z ) -> (s y , -s z ) (16) 

t , (S X , Sy , >S 2 ) > ( S X , Sy , S Z ^j (17) 

~t , Sy -> Sy (18) 

T 

t + — , (s x , S g ) -> (-Sa, , -S z ) (19) 

If we add dissipation terms, these terms will break time reversal symmetry, and we are 
left only with ([TB] ) and flT^) . All of the above statements for the quantum system can be 
recovered. Especially nonzero dissipation and 7, a 7^ lead to nonvanishing magnetization 
along the y-axis, even for MCF. 

Let us summarize the presented results. We have shown that time-periodic magnetic 
fields with zero mean may induce nonzero averages of spin components which would be 
strictly zero in the absence of these fields. The spin is simultaneously experiencing some 
local anistropy field or simply an external constant field. In addition the spin is coupled 
to some thermal environment characterized by some finite temperature and a characteris- 
tic relaxation time ||. The reasoning follows symmetry considerations of the dynamical 
equations. In the case of a classical spin these equations formally coincide with the Heisen- 
berg equations for the quantum spin operators. In the quantum case we instead solve the 
(purely linear!) equations of motion for the independent components of the density matrix. 
Remarkably the symmetry properties obtained from both approaches coincide. 

The quantum approach shows that for infinite temperatures all spin component averages 
will vanish. This follows from p(t — > 00; f3 — > 0) = and Tr S Xi y jZ = 0. 

For the spin 1/2 case we proposed a MR experiment to observe the effect. One should 
choose the time-periodic magnetic field to be not perpendicular to the static magnetic field. 
Further the amplitude of the time-periodic field should be not too small such that the 
generated S x and S y components are measurable. The attenuation spectrum should show 
resonances located at subharmonics of the original resonance. The intensity of the sattelite 
peaks is a function of both the angle between both fields and the inverse relaxation time v. 
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Experiments which probe directly the nonzero spin components can be performed by 
adding yet another probing field to the system, and varying its frequency while keeping the 
frequency of the original probing field. This will be studied in detail in future work. 
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Figure captions. 

Fig.l. 10S X (solid), lOSy (dashed) and S z (dotted) as functions of cu (see text for pa- 
rameters) . 

Inset: S XjVjZ versus time for one period of h(t) at u — 1.5 (same line codes as in Fig.l). Note 
that functions are not scaled here! 

Fig.2. 

S y as a function of uj (see text for parameters). 
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